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Abstract 



Using the two-dimensional Jordan- Wigner fermionization we calculate the thermodynamic quan- 
tities of the (spatially anisotropic) square-lattice spin-| anisotropic XY (XZ) model. We compare 
^ the results of different approaches for the ground-state and thermodynamic properties of the model. 
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Two-dimensional (2D) quantum spin models have been extensively studied during last years mainly 
because it is believed that they may be of use for describing the magnetic properties of Cu02 layers in 
the high-temperature superconductors JjJ]. There exist a number of analytical approaches for a study 
of the thermodynamic properties of 2D quantum spin models, e.g., the conventional spin-wave analysis, 
the Green function technique, the approach based on the 2D Jordan- Wigner fermionization as well as 
the coupled cluster method, the correlated basis function method etc.. In what follows we consider the 
spin-i anisotropic XY model on a spatially anisotropic square lattice within the framework of the scheme 
based on the 2D Jordan- Wigner fermionization and compare the results derived for the ground-state and 
thermodynamic quantities with the exact ones (ID limit, square-lattice Ising model) and the predictions 
of other approximate theories. The performed calculations yield an impression about the region of validity 
of some approaches usually applied for a study of thermodynamics of 2D quantum spin models. 

We start from a model of N — > oo spins \ on a spatially anisotropic square lattice governed by the 
anisotropic XY Hamiltonian 



H = E E J ( (1 + y) s tj s t+i,j + a - 7)<^r+i, 

i=0 j=o \ 

+Jx((l +7)<i*? J - + i + (1 - -<.,<,..,))• 



(1) 
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Here J and J± — RJ are the exchange interactions between the neighbouring sites in a row and a column, 
respectively (for concreteness we assume both to be positive) , and the parameter 7 controls the anisotropy 
of the exchange interaction. Making use of the 2D Jordan- Wigner fermionization and adopting a mean- 
field treatment of the phase factors which appear after the fermionization j^, |3j we perform consequently 
the Fourier and Bogolyubov transformations to arrive at the following Hamiltonian of nonintcracting 
spinless fermions which represent the initial spin model ([!]) 



# = ELE A «( k ) (<«v<,«-~) 



(2) 



Ai(k) = sj (Jj_ cos k y + ~fJcosk x ) + (Jsirik x + jJ± sinfc. 



VJ > 



A2(k) = y (Jj_ cosfc a — 7 J cosfc x ) + (Jsinfc x — 7J1 sinfcj,) 

(the prime denotes that k in the thermodynamic limit varies in the region —n < k x < 7r, —it + \k x \ < 
k y < 7r — \k x \). The Helmholtz free energy per site 



f= ~wL^ L ^ l ln v cosh 2 ) + ln I 2 cosh ^ J J (3) 

yields the thermodynamic properties of the spin model (Q). In Fig. 1 we plot the ground-state energy 
per site of the spin model (Q), (0) (dotted curves) in comparison with the exact results if R — (ID 
model) or 7 = 1 (square-lattice Ising model) and the spin- wave theory result for 7 = 0, R = 1 (spatially 
isotropic square-lattice isotropic XY model). Eq. @ contains the exact result in ID limit (Fig. lb), 
however, deviates noticeably from the exact result for 7 = 1 (compare the curves 3 in Fig. la). For 7 = 0, 
R — 1 Eq. (|§|) yields the result which differs from the spin-wave theory prediction denoted by the full 
circles. (The outcomes of different numerical approaches (see Q]) lie within the full circles.) From the 
exact calculation for 7 = 1 J5j we know that the temperature dependence of the specific heat exhibits a 
logarithmic singularity. Obviously, the Jordan- Wigner fermions (^|), (|3|) cannot reproduce this peculiarity 
inherent in the spin model. 

It is worth to remind here that the conventional spin-wave theory was originally thought to be un- 
satisfactory for quantum XY models However, the authors of the paper |Q showed that considering 
the XZ rather than the XY Hamiltonian one gets within the spin-wave theory satisfactory results of the 
same quality as for the Heisenberg Hamiltonian. Following this idea we perform the rotation of the spin 
axes s x — > — s z , s y — > s x , s z —s v and consider instead of ([!]) the following Hamiltonian 

00 00 j 

H = E E ( J { {1 ^ s Zi s i+i,j + a + 7)<i^u 

i=0 j=0 ^ 



+■71 ((1 - iKAj+i + (1 + 7)<i<i+i) ) • ( 4 ) 



Proceeding further with (|J) in the described above manner and assuming (for concreteness) antiferromag- 
netic long-range order while decoupling the quartic fermionic terms || we get instead of (||) the following 
Hamiltonian 

' 2 / 1\ 

F = ElE A «( k ) )+N(l + ^)(J+J±)m 2 , (5) 

a=l ^ ' 



Ai, 2 (k) = 2yU 2 + £ 2 + C 2 + P 2 + M 2 ± lyJlABCV + ^ 2 P 2 + B 2 C 2 + M 2 (B 2 + £> 2 ), 

1 — 7 1 — 7 

A = J±_ cos k v , B = J± sin k v , 

4 4 

C = - — — J sin k x , T> = - — — J cos k x , 

4 4 x ' 

M = (1 + 7 )(J + Jj_)m, 



2 



where m is determined self-consistently by minimizing the Hclmholtz free energy per site 

^J- 7! J- K {2^f\ dm 2 9m 2 J 

In Fig. 2 we plot the ground-state energy of the spin model (|J), (|5|), (^) (dashed curves). The results 
based on Eqs. (g), (||) for 7 = 1 reproduce the exact result for square-lattice Ising model (curve 3 in Fig. 
2a) as well as the spin-wave theory prediction for 7 = 0, R = 1. However, the result based on Eqs. (||), 
(H) for i? = does not coincide with the exact one in ID limit (curve 1 in Fig. 2b). 

To summarize, we have calculated the thermodynamic quantities for the spin-^ anisotropic XY (XZ) 
model on a spatially anisotropic square lattice using the 2D Jordan- Wigner fermionization. To reveal a 
quality of the results obtained within the framework of this approach we have compared them with the 
exact results available in ID limit and extremely anisotropic exchange interaction limit (Ising interaction). 
We have found that although there is an agreement with the spin-wave theory and other approximate 
approaches a disagreement with the exact results may be noticeable. Thus, the question about the quality 
of the results based on the 2D Jordan- Wigner fermionization (as well as, e.g., of the spin- wave theory 
results) remains still open and requires further studies. Moreover, for arbitrary values of anisotropy 
parameter 7 (7 ^ 1) for R ^ a comparison with the exact diagonalization data and results of other 
numerical approaches is desirable. 
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partly supported by the DFG project 436 UKR 17/101. O. D. acknowledges the kind hospitality of the 
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FIGURE CAPTURES 

FIGURE 1. The ground-state energy per site for the square-lattice spin-i anisotropic XY model (Q) 
e vs. R (a) (1 - 7 = 0, 2 - 7 = 0.5, 3 - 7 = 1) and e vs. 7 (b) (1 - R = 0, 2 - R = 0.5, 3 - R = 1); 
exact results (solid curves) and the approximate results obtained on the basis of (Q) (dotted curves); the 
full circles correspond to the spin- wave result for 7 = 0, R = 1. 

FIGURE 2. The same as in Fig. 1 for the XZ Hamiltonian (Q); the approximate results obtained on 
the basis of (|J), (^|) are shown by dashed curves. 
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